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By H. F. Poppendiek 

ABSTRACT 

An analytical heat transfer solution is derived and evaluated for 
the general case of a turbulently flowing liquid metal which suddenly 
encounters a step -function boundary temperature in a channel system. 
Local Nusselt moduli, dimensionless mixed-mean fluid temperatures, 
and arithmetic -mean Nusselt moduli are given as functions of Reynolds 
and Prandtl moduli and a dimensionless axial -distance modulus. These 
solutions are compared with known solutions of more specific systems 
as well as with a set of experimental liquid-metal heat transfer data for 
a thermal entrance region. 

^Originally prepared as Report ZPh-015, by the Physics Section of 
CONVAIR, a division of General Dynamics Corporation. Reproduced in 
the original form by NASA, with CONVAIR permission, to increase 


availability. 
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NOMENCLATURE 


SYMBOLS 


a, 



f(r) , 


am 


h 


x* 


J o’ Y 


o’ 


J’ ,Y' , 
o o 

k. 


P. 


2 

thermal molecular dlffusivity, k/yc^, ft /hr 

series coefficients in equation (14) defined by equation (15), 
dimensionless 

half the channel breadth (see Figure 1) , ft 

series coefficients in equations (11) and (13), dimensionless 

constant in equation (37) 

constant in equation (30) , dimensionless 

constants in the general solution, equation (9), dimensionless 

heat capacity, Btu/lb°F 

a function of the variable r, dimensionless 

arithmetic mean unit thermal conductance or heat transfer 

coefficient, -9-/A t , Btu/hr ft 2o F 
A am’ 

local unit thermal conductance or heat transfer coefficient, 

(k)/<‘w - l m >- Btu/hr f ‘ 2 ” F 

zero order Bessel functions of the first and second kind, 
respectively 

derivatives of the zero order Bessel functions of the first 
and second kind, respectively, with respect to the argument 

thermal conductivity, Btu/hr ft^ (°F/ft) 

constant in equation (37) 


S, 

t. 


t . 
m 


2 

local heat flux, Btu/hr ft 

a function of x and y given by equation (37), dimensionless 
liquid-metal temperature at some point x, y, °F 
mixed-mean fluid temperature, "F 
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SYMBOLS (Continued) 


uniform wall temperature (step function), °F 
initial liquid-metal temperature (see Figure 1), °F 

liquid-metal velocity profile (function of y) , ft/hr 

mean or uniform liquid-metal velocity, ft/hr 

sum of J and Y functions in equation (11) 
o o 

independent variable defined in equation (39) , dimensionless 


Cartesian coordinates, ft 


(t - 1 ) + (t - t ) 


arithmetic-mean temperature difference, , °F 

2 

2 

eddy diffusivity, ft /hr 

2 

kinematic viscosity, ft /hr 

parameter arising in the separation -of -variables technique, dimensionless 
eigenvalues of equation 12, dimensionless 

3 

weight density, lb/ft 
gamma function 


DIMENSIONLESS MODULI 

4 

F = 
o Pe 


F = 

1 0.1 

Re 


h 4b 
„ X 

Nu^ - local Nusselt modulus 

h 4b 
am 

Nu = , arithmetic -mean Nusselt modulus 

am k 

Pe = Re . Pr , Peclet modulus 


yc v 

p 

Pr = , Prandtl modulus 

k 




DIMENSIONLESS MODULI (Continued) 

U 4b 

Re = — , channel Reynolds modulus 
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INTRODUCTION 

It is well known that the convective heat transfer in the entrance regions of duct 
systems where thermal and hydrodynamic boundary layers are not yet established can 
be far superior to heat transfer in the established flow regions. A quantitative under- 
standing of this type of heat transfer, sometimes called entrance region heat transfer, 
is essential when designing high heat-flux cooling systems for rocket motors, nuclear 
reactors, exhaust nozzles, and missile nose cones. Because the liquid metals are the 
most effective high -temperature coolants known, they are considered exclusively in 
this report. Further, in many practical flow systems the hydrodynamic boundary layers 
have been completely or almost completely established before the thermal entrance 
region is encountered. Therefore, the work presented deals with heat transfer in the 
thermal entrance region with an established velocity field. 

A thermal entrance region results when a thermally established fluid, flowing in a 
duct of uniform cross section, suddenly encounters duct surface with some new boundary 
temperature distribution. Under these circumstances the temperature field in the fluid 
is no longer established and thus greatly influences the local convective heat transfer. 

As an example, for a step function boundary temperature entrance region the heat 
transfer can be very high because the thin thermal boundary layers have low thermal 


resistances. 
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A limited number of turbulent flow entrance region solutions for ducts are avail- 
able in the literature. Sanders (1) has obtained a turoulent flow solution for a step 
function entrance region in a pipe by transforming th3 turbulent core to a laminar 
core of equivalent thermal resistance. Seban and Shimazaki (2) have obtained some 
specific numerical solutions for the uniform wall -heat-flux entrance region in a 
pipe. Elser (3) obtained a simplified solution for only the initial portion of the ther- 
mal entrance region of a pipe containing a turbulent Tuid and a step function wall 
temperature distribution. Several mathematical ana.yses for forced convection heat 
transfer in thermal entrance regions for low Prandtl modulus (liquid metal) systems 
were presented previously (4); three dealt with low, turbulent Reynolds moduli (radial 
heat flow by conduction only) and two others dealt with all turbulent Reynolds moduli 
(radial heat flow by eddy transfer as well as conduction). The solutions for the latter 
two problems had been completed but not evaluated. One solution has since been evalu- 
ated and the results are presented. The solution desc ribes the general case of turbu- 
lently flowing liquid metal which suddenly encounters a step function boundary temperature 
in a channel system. The derivation of the mathemat cal solution is presented. Local 
Nusselt moduli, dimensionless mixed-mean fluid temperatures, and arithmetic-mean 
Nusselt moduli are given as functions of Reynolds and Prandtl moduli and a dimension- 
less axial-distance modulus. The solution is shown to reduce correctly to known specific 
solutions of the general case. Also the solution is compared with a set of experimental 
liquid-metal heat transfer data previously obtained in a thermal entrance region. 
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ANALYSIS 


The idealized system which defines heat transfer in a liquid metal flowing 
turbulently through a channel (between two parallel plates of Infinite extent) with a 
step function boundary temperature is based on the following postulates: 

1) The wall temperature distribution is a simple step function, 

t = t G for x < 0 and t = for x >0. The fluid approaching the en- 
trance region has an established temperature, t Q (see Figure 1). 

2) Longitudinal heat conduction is small compared to convection and 
is neglected (see Appendix 1). 

3) The established turbulent velocity profile is represented by a uni- 
form distribution u = U, (see Appendix 2). 


4) 


5) 

6 ) 


The eddy diffusivity distribution varies linearly with distance from 
the wall and as the nine tenths power of the Reynolds modulus, 


= C± Re^ ' 9 i , (see Appendix 3) . 
b 


The fluid properties are invariant with temperature. 
Steady state exists. 


The differential equation describing the convective heat transfer in dimension- 
less form is (see Appendix 3) 


d T 

ax 



ay 


(i) 


This equation can be expressed in a simpler form by making the change of variable, 


r = \ ~ 4 Y = a|f + F Y' 

y Pc o.i ' ° 1 

* Re 


( 2 ) 
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The resulting boundary value problem to be solved is, 


F 2 2 F 2 

1 9 T 1 1 9T 

(3) 

4 a r 2 4 r 9r 

T (o, r) = 1 

(4) 

TiX.-fF) =0 

(5) 

ar' ( MV V-° 

(6) 

lim T (X, r) = 0 
X- <*> 

(7) 


This problem can be solved by the separation-of-varlables technique. Let 


T = 


^(X), * 2 <r) 


( 8 ) 


where ^(X) and 4> 2 (r) are functions of X and r, respectively. Upon substituting 
equation (8) into equation (3), two total differential equations result. One involves 
♦jCX) and the other 4> 2 (r). Their solution yields, 


1 o \F^ / 2 


(9) 


where p is the parameter arising in the separation -of-variables technique and c 


and c are constants in the general solution. 

Z 


From equation (5), the constant, 


V 


is found to be 
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(10) 



The constant p can be evaluated by substituting the temperature function, T, 
Into the boundary condition given by equation (6). The resulting expression is, 



which is the eigenfunction. The terms P are the eigenvalues (n = 1, 2, 3 . . . ) 

The constant c^ in equation (11) is now replaced hy the constants b (n = 1, 2, 3, 
corresponding to the values, P ^ . The constants, b , can be evaluated from the boundary 
condition given by equation (4), 


-E 

n = l 



(13) 


It can be shown (see Appendix 4) that a function f(r) can be expanded Into a U 

o 

series over the interval r tor , 

1 2 
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00 

' <r> ■ £ \ u oGf ) 

n = 1 1 


where 



( 14 ) 


(15) 


For the boundary value problem being considered f(r) = 1. From equations (13) and 
(14), 



Thus, the temperature solution Is, 



where the coefficients A^ are given by equation (15) and the eigenvalues are defined 
by equation (12). 

As X approaches zero, many terms are required to obtain convergence in the 
series solution given by equation (17). Thus, It Is convenient to ubo an asymptotic 
solution In that region (see Appendix 5). 
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The mixed-mean fluid temperature , t , car be expressed in dimensionless 

m 


form as 


Vf + F ’ 
X O 1 


T m ' Tfdr 

K 


The local Nusselt modulus, Nu^ # can be expressed as, 


x k 


<3Y/Y~0 




where 


00 -p X 28 




* j oVr 'i v Y 


oo .2 

-P X 2P / F 

V -1 a _ « _JL I 2 1 


IrffJCv - ^ 


3 A pr 

ny O 


Upon substituting equations (20) and(18) into equation (19), there results 


1. The mixed-mean fluid temperature for a channel system is defined as 

r b 

/ u t dy 

Jo 


l m <~b 
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The arithmetic mean Nusselt modulus, Nu , which Is based on an arithmetic 

am 

mean wall -fluid temperature difference , is expressed as 


Nu 

am 


h 4b 
am 

k 


( 22 ) 


where 


h 

am 


= a 

A At 

am 


At 

am 


(t - t ) + (t - t ) 
wo w m 


(23) 


(24) 



(25) 


Upon substituting equations (23), (24), and (25) into equation (22), there results 


Nu 


am 



dX 


2 


(26) 



RESULTS 


The eigenvalues and series coefficients for this boundary value problem 

were evaluated for a wide range of Reynolds and Prandtl moduli. The quantities 

T , Nu y , and Nu were calculated as functions of the parameter X, and Rey- 
in am 

nolds and Prandtl moduli. The results are presented in Figures 2 to 20 and in 
Appendix 6 (Tables I to VI). From the Nu^ graphs it may be observed that the 
entrance length increases as Reynolds modulus increases. The established 
Nusselt moduli, Nu , are shown as a function of Reynolds and Prandtl moduli, in 

QO 

Figure 8. The T m graphs show that the fluid temperature approaches the wall 
temperature in a shorter distance from the entrance as the Reynolds modulus de- 
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Figure 5 .- Nu y vs X and Re for Pr = 0.02. 
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Figure 7-- Nu x vs X and Re for Pr = 0.06. 
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Figure 9.- T m vs X and Re for Pr = 0.002 



23 



Figure 10.- T m vs X and Re for Pr = 0.005 
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Figure 11.- T m vs X and Re for Pr = 0.01. 
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Figure 13- - T m vs X and Re for Pr = 0.03 
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Figure 14.- T m vs X and Re for Pr = 0.06 
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Figure 15 . - Nu am vs X and Re for Pr = 0.002. 
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Figure 16.- Nu^- vs X and Re for Pr = 0.005- 
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Figure 17- - Nu arn vs X and Re for Pr - 0.01 
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Figure 18.- Nu™ vs X and Re for Pr = 0.02. 
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DISCUSSION 


Comparison of Solution with Specific Analytical and Numerical Work 

The well known solution for the special case of flow in a channel under low Rey- 
nolds and Prandtl moduli conditions (radial heat flow by conduction only) with a uniform 
velocity profile was compared with the general convection solution derived above for a 
condition where eddy transfer is small compared to conduction, namely, Re = 5000 and 
Pr = 0. 002; the respective temperature and Nusselt modulus solutions for the two cases 
were in complete agreement for all values of X. 

The established Nusselt modulus solution for long channels, Nu^ for uniform 
wall -heat-flux conditions by Martinelli (5) was also compared to the established Nusselt 
modulus results (Figure 8) obtained from the general solution* Figure 21 illustrates that 
the two sets of calculations are in good agreement, anc the uniform wall temperature 

2 

Nusselt moduli fall about 15 percent below the uniform wall -heat-flux Nusselt moduli. 

A discussion of the linear approximation of the complicated eddy diffusivity pro- 
file was presented in Appendix 3. It was concluded tha: the agreement between the linear 
and actual eddy diffusivity functions in the important heat transfer layers was good. To 
substantiate this conclusion further, a numerical analysis of the convective heat transfer 
problem was made with the actual eddy diffusivity function for the specific case of Re = 

200, 000 and Pr = 0.01. The agreement between the analytical solutions (using the linear 
eddy diffusivity function) and the numerical solutions (tsing the actual eddy diffusivity 
function) was good. For example, the local Nusselt moduli for the two cases shown in Fig- 
ure ZZ fall within about 10 percent of each other. 

2. It has been found (6) that for the case of established flow X -+~<*>the Nusselt moduli for 
uniform wall-heat-flux boundary conditions are a little greater than those quantities for 
uniform wall -temperature boundary conditions; this result was based on analytical conduc- 
tion solutions as well as on numerical turbulence solut ons. 



Figure 21.- Nu^ vs Re and Pr for uniform wall heat flux and uniform wall temperature cases 
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Com] ) ; irison of Solution with Experimental Data 

The only experimental heat transfer data that could be found for a uniform wall- 

temperature system were those of Harrison (8) for turbulently flowing mercury in small 

diameter pipes with short heating sections. These experimental data occurred in the 

3 

initial part of the entrance region where the asymptotic solutions (Appendix 5) hold. 
These solutions compared favorably with Harrison’s experimental results (4) as shown 
in Figure 23. 

Application 

There are a number of practical high-heat-flux cooling systems to which the con- 
vection solution presented may be applied. One good example is a nuclear reactor core 
whose fuel elements are plates which are cooled by a liquid metal. It is possible to cal- 
culate ideal fission heat source distributributions in the axial direction in fuel plates of 
nuclear reactors which yield maximum uniform fuel plate temperature distributions. 
Maximum reactor powers would result from such axial temperature profiles. 

Another example is a liquid-metal cooled target of an accelerator. For high-con- 
ductivity targets, the power that can be absorbed without exceeding limiting target tem- 
peratures can also be determined. 

The solution is also applicable to special cooling problems encountered with missile 
nose cones and exhaust nozzles of propulsion systems; the excellent heat transfer that can 
be obtained in the thermal entrance regions of liquid-metal convection systems can be 
used to advantage. In such cases, the liquid metal would flow between plates (fins) which 
would be attached to the heat transfer surface to be cooled. 

3. Mean Nusselt moduli, Nu am , in Harrison’s system were from about two to four times 
higher than those for established flow values, Nu^. 



EXPERIMENTAL Ni 



Figure 23.- A comparison of Harrison's exper:.mental mercury heat transfer 
data with the asymptotic solution for a p:.pe with a uniform velocity 
profile . 
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APPENDIX 1 

LONGITUDINAL HEAT CONDUCTION 

Mathematical convection solutions are usually based on the postulate 
that the longitudinal heat conduction term is negligible compared with the 
convection term. The validity of this premise was investigated for channel 
and pipe systems for low Prandtl modulus conditions and uniform velocity 
profiles (reference 4). Temperature solutions were derived for the case 
where the axial conduction term was not neglected and the results were 
compared to the solution for the case where axial conduction was neglected. 
The comparison revealed that axial conduction was relatively unimportant 
for turbulent flow conditions. For example, it was found that for a Prandtl 
modulus of 0. 005 (which represents a practical minimum value for liquid 
metals) and the low turbulent Reynolds moduli s of 8,000, the two Nusselt 
moduli for a pipe system differed by 2. 5 percent at 0.4 of a diameter from 
the entrance, 1. 3 percent at one diameter from the entrance, and a still 
smaller percentage at greater distances from the entrance. 
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APPENDIX 2 

FLUID VELOCITY DISTRIBUTION 


The blunt-nosed turbulent velocity distribution in a channel can be 
represented satisfactorily by a uniform distribution. In reference 4, entrance 
region heat transfer solutions for a channel were presented for low turbulent 
Reynolds modulus conditions for both uniform and blunt-nosed (one seventh 
power law) velocity distributions. The results are graphed in Figure 24. 

4 

Note that the Nusselt moduli for these two profiles differ from each other 
by percentages varying from about 6 to 20 over the wide range of Re Pr/X 
values shown. The higher difference-percentages correspond to the initial 
portion of the entrance region which normally represents only a small 
fraction of the total heat transfer surface in the entrance region. However, 
if one desires to calculate convective heat transfer in the very initial portion 
of the entrance region, where a non-uniform velocity profile should be used, 
the asymptotic solutions given in reference 4 can be used. At the higher tur- 
bulent Reynolds moduli, the actual turbulent velocity profiles become very 
flat (reference 7); in this region the idealized uniform velocity profile very 
closely represents the actual ones. 


4. The results for a third velocity profile, the parabola, are also graphed 
in Figure 24 for purposes of comparison. 




Figure 24.- Comparison of low Reynolds number (conduction) heat transfer solutions for uniform 

one seventh power law, and parabolic velocity profiles. 
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APPENDIX 3 

EDDY DIFFUSIVITY DISTRIBUTION 


The general differential equation describing convective heat trans- 
fer for the idealized system described previously Is 


U 


jH 

9x 


9 y 


(a + €) 


at 

a y 


This equation can be expressed in a dimensionless form. 


9 T 9 4/1 | (\ >T 

9X 9 Y Re \J>r + v) 9Y 


(27) 


(28) 


The analogy between heat and momentum transfer has been firmly established 
in many experimental turbulent-flow systems. Thus it is postulated in this 
analysis that heat transfer and momentum transfer eddy dlffuslvities are 
identical . Momentum transfer eddy diffusivlties (in dimensionless form) 
for a channel system can be represented as follows: 


Laminar Sublayer 

, 0. 9 

(0 < Y < 131.5/Re ) 

Buffer Layer 

, 0.9 , 0.9 

(131. 5/ Re < Y < 789/ Re ) 

Outer Turbulent Layer 
0.9 

(789/ Re < Y <0.5) 

Inner Turbulent Layer 
(0.5 < Y < 1.0) 


€ 

V 

€ 

V 

C 

V 

€ 

V 


Y 


0. 9 y 

0.0076 Re ~ - 1 
b 


X29) 


0.9 y v 
0.0152 Re (1--)- 
b b 


0.0038 Re 


0. 9 


J 


5. These relations are obtained from the shear stress expression, the 
generalized velocity profile for turbulent flow, and the shear equation. The 

results pertain to smooth channels over a Reynolds modulus range of 

3 6 

5xl0 to 10 . 



b6 


Because the important heat transfer layers are those nearest the wall, an 
idealized eddy diffusivity function is postulated which approximates the 
actual one in that region; the idealized eddy diffusivity relation is, 

€ . 9 y 0 . 9 y 

— = C Re -f" = 0.01 Re -f- (30) 

v 1 b b 

Note, from equation (28), the dimensionless eddy diffusivity always appears 
together with the reciprocal of the Prandtl modulus in the form of a sum, 

. A comparison of this sum for the actual and idealized eddy dif- 

Pr v 

fusivity functions is shown in Figure 25 for a typical liquid- metal system. 

Note, that in the important outer half of the flow channel, the idealized sum, 
1 « 

— + — , is a good approximation of the actual quantity. Upon substituting 
Pr v 

equation (30) into equation (28), equation (1) can be obtained. 


6. The heat transfer layers between the wall and about half the distance 
to the duct center are the important ones because a large fraction of the total 
radial temperature drop is found there. This is true because 1) the radial 
heat flows and 2) the thermal resistances are large in this region in compar- 
ison to the central duct core. 
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APPENDIX 4 

SERIES EXPANSION OF f(r) 


U 


In equation (13) it was found necessary to expand unity into series of 
— r) functions. The known procedure for doing this (see reference 


9 for example) is discussed here for the case of a general function, f (r). It 

/ 2 P n \ 

is desired to expand f(r) into a U — r) series in the interval r tor , 


°V F 1 / 

oo /2 p 

f (r) = Y' A U ( — 11 

n = l n °\ F l 


(31) 


"2p 


If both sides of equation (31) are multiplied by the term, r U ( — — r) 

\ F 1 J 

and the equation is integrated from r to r , all integrals involving 

X 2 

/ 2 P n \ 

U l r ) terms with two different values of p are zero. Only the 
X / /2P n v 

integral having U ( ~ r) terms with the same values of p are not 

°\ F i / n 

zero. Thus the equation for the series coefficient is found to be 

r 



(32) 


The proof that the integral 


U 


'2p 1 
n 


dr 


is equal to zero is outlined below. The quantities p and p* represent two 

n n 

different eigenvalues. The proof usually involves writing down two differential 

Bessel equations for two different solutions, one in terms of p and another 

n 

in terms of p^. Upon multiplying each equation by the solution of the other, 
subtracting, and integrating over the range r to r , there results 



4 9 



Upon substitution of equations (5) and (12) into the bracket of equation (33), the 
integral is found to be, 






50 


where the limits for the specific problem being studied are 

r = \[f + F and r = Vf^ 

2 * o 1 1 o 
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APPENDIX 5 

ASYMPTOTIC CONVECTION SOLUTION 


Equation (17) converges very slowly for small values of X. A general 
asymptotic solution has been derived (reference 4) which can be used to 
evaluate the temperature and heat transfer in this region. For small values 
of X, the turbulent flow system being studied here reduces to the case of 
convection over a single flat plate with radial heat flow being achieved entirely 
by conduction. This is true because for small values of X 1) the influence of 
heat flow from the other wall would not exist and 2) the thermal boundary 
layer would not have diffused into the turbulent flow region. The boundary 
conditions for the asymptotic case are 


lim 

X*«o 


lim 
Y - « 


T (X, 0) 
T (X, Y) 
T (0, Y) 
T (X, Y) 


0 

0 

1 

1 


A 


J 


(36) 


The solution of equation (1) together with the boundary equations (36) can be 
accomplished by making a change in variable. 


S 



(37) 


where p and c are constants to be determined. Upon substituting equation 
(37) into equation (1), the solution for the boundary equations (36) is found 
to be (reference 10) 



(38) 
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where 


w 


Ub 
4 a X 


(39) 


Further It can be shown that for small values of X, 


and 






Nu 

am 


2 Re 


1/2 1/2 

Pr 




(40) 


(41) 


CONVAIR, 

A Division of General Dynamics Corporation, 
San Diego, Calif., May 1958- 
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